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Abstract We discuss the HFB equations in coordinate representation, a suitable 
method for handling the full effects of the continuous quasiparticle spec- 
trum. We show how the continuum HFB equations can be solved with 
the correct asymptotic conditions instead of the discretization condi- 
tions which are commonly used in the literature. The continuum HFB 
method is illustrated with a model where the mean field and pairing 
field have simple forms. The relationship with the continuum Hartree- 
Fock-BCS (HF-BCS) approximation is also discussed. Realistic HFB 
and HF-BCS calculations based on Skyrme interactions are compared 
for the case of a neutron-rich nucleus. 



1. INTRODUCTION 

In many nuclear systems pairing correlations have an important in- 
fluence on physical properties. This situation occurs when the energy 
difference between an occupied orbital and a neighbouring unoccupied 
one is relatively small, thus enabling a nucleon pair to be promoted to 
the unoccupied level by the interaction among the pair. Well known the- 
oretical approaches such as the Hartree-Fock-Bogoliubov (HFB) or the 
Hartree-Fock-BCS (HF-BCS) approximations have been developed [1] 
to treat the pairing effects in nuclei. 

In usual situations, i.e., when the nuclear system is stable and far 
away from the drip lines the active orbitals are well bound and it may 
be sufficient to solve the HFB or HF-BCS equations within the discrete 
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subspace of those active orbitals. When one approaches a drip line this is 
no longer true since the active orbitals must also include states belonging 
to the continuous single-particle spectrum. One of the most convenient 
tools for the treatment of pairing correlations in the presence of contin- 
uum coupling is given by the HFB approach. The standard procedure is 
to discretize the quasiparticle continuum by solving the HFB equations 
in a finite basis of orthonormal functions, or by imposing to the solutions 
of HFB equations, written in coordinate space, a box boundary condi- 
tion at some distance i?& [2, 3, 4]. There may be circumstances where 
a discretized solution is unsatisfactory unless an extremely large box ra- 
dius Rb is used. Such should be the case for halo nuclei, for instance. 
In these lectures, we show how the HFB equations can be solved with 
correct boundary conditions. Based on these solutions we examine how 
much the continuum, especially the resonant continuum, could affect the 
pairing properties of nuclei close to the drip lines. 

The structure of the lectures is the following: In the next section we 
introduce the HFB equations in coordinate space [5, 2] and in section 3 
we examine the solutions which satisfy the correct boundary conditions, 

1. e., of scattering wave type. In section 4 these considerations are illus- 
trated for the particular case of a square well potential and a constant 
pairing field of a finite range. Then we discuss how the resonant contin- 
uum can be taken into account in HF-BCS approach. Finally we show 
how the continuum, calculated in different approximations, can affect 
the pairing properties of a nucleus close to the neutron drip line. 

2. HARTREE-FOCK-BOGOLIUBOV 
EQUATIONS IN COORDINATE SPACE 

Let us denote by H the hamiltonian of the nuclear system. When 
treated in the HFB approximation, H will give rise to the mean field 
and the pairing field in which the quasiparticles are moving. We assume 
that the two-body effective interactions associated with the Hartree-Fock 
(HF) mean field and pairing field are zero-range forces so that the total 
energy E is a functional of the local particle density p(r) and pairing 
density rc(r): 

p(r) = {HFB\V + (r)V(r)\HFB) , 

«(r) = {HFB\^ + (r)^ + (r)\HFB) , (1.1) 

where v I' + (r) is a nucleon creation operator and \HFB) is the HFB 
ground state. Typically, the self-consistent mean field is generated by 
a Skyrme-type interaction whereas the pairing field is produced by a 
zero-range, possibly density-dependent interaction. 
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The HFB equations in coordinate space can be expressed in the fol- 
lowing form: 

(h-X A \ (U a (r) \_ (U a (r) \ 

[ A -(h-X) )[v a (r) )- Ea {v a (r) ) ^ 

where A is the chemical potential, h is the HF hamiltonian and A is 
the pairing field. The local nucleon density and pairing density can be 
written in terms of the solutions (U a ,V a ) in the form: 

p(r) = J2\V a (r)\\ (1.3) 



«(r) = J>«(r)l£(r). (1.4) 

a 

For Skyrme forces, the HF and pairing fields may also depend on deriva- 
tives of p and k and on some local currents because of the velocity de- 
pendence of the interaction. Since all densities appearing in the energy 
functional are local, Eq.(1.2) is a set of coupled differential equations 
which is highly non-linear because the fields h and A depend themselves 
on the solutions U a ,V a (self-consistency problem). Nevertheless, using 
Skyrme-type forces brings a major simplification. If one starts instead 
with a finite range effective force, for instance a Gogny's force, Eq.(1.2) 
would be a set of coupled integro-differential equations which is usually 
solved in a harmonic oscillator basis. The general properties of Skyrme- 
HF hamiltonians are well-known and detailed expressions of h and A for 
the case of Skyrme forces can be found for instance in Ref. [2]. 
The Skyrme-HF hamiltonian has the general form: 

h = - v -d^) v + ^ (r) - <L5) 

where the effective mass m* and HF potential Vhf depend on nucleon 
densities and currents. An important property is that m* tends to the 
nucleon mass m and Vh f tends to zero at infinity, with the same rate as 
the densities go to zero. We shall use this property in the next section 
to establish the asymptotic forms of the solutions of Eq.(1.2). 

For illustration, we can look at the pairing field calculated with a 
pairing interaction often used in the literature [3, 4]: 



V(l,2) = V 
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<S(ri - r 2 ) . (1.6) 
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To be meaningful, such an interaction must be used with an energy cut- 
off in the quasiparticle spectrum, i.e., the summations in Eqs. (1.3,1.4) 
must be limited to E a < E cuto ff- In this case, the pairing field is: 



On this example we can see that the local pairing field behaves asymp- 
totically like the local pairing density. 



3. ASYMPTOTIC BEHAVIOUR OF HFB 
SOLUTIONS 



From the symmetries of Eq.(1.2) it can be seen that, if (E a ,U a ,V a ) 
is a solution then (— E a , V*, U*) is also a solution. Thus, we need to 
consider only one class of solutions and we choose those with E a > 0. 

The asymptotic behaviour of the solutions of Eq.(1.2) has been dis- 
cussed in detail by Bulgac [5]. Let us consider for simplicity the case of 
spherical symmetry and write: 



Eq.(1.2) now becomes a one-dimensional equation in the radial coordi- 
nate r for each a = At very large distances the HF hamiltonian 
h tends to —(h 2 /2m)(^-^jr — ) (plus a Coulomb potential for pro- 
tons) while the pairing field A(r) has vanished. The equations for u a (r) 
and v a (r) are decoupled and one can easily see how the physical solu- 
tions must behave at infinity. Thus, for a negative chemical potential 
A, i.e., for a bound system, there are two well separated regions in the 
quasiparticle spectrum: 

- between and —A the quasiparticle spectrum is discrete and both 
upper and lower components {u a (r),v a (r)) of the HFB wave function 
decay exponentially at infinity; 

- above —A the quasiparticle spectrum is continuous and the physical so- 
lutions are such that at infinity the upper component of the HFB wave 
function has a scattering wave form (see the next section) while the 
lower component is exponentially decaying. In what follows the contin- 
uous HFB wave functions are normalised to a delta function of energy. 

Thus, the summations in Eqs. (1.3-1.4) should in fact include integra- 
tions over the continuum of the quasiparticle spectrum: 



A(r) 




(1.7) 




(1.8) 



P {T) 



0<E a <-\ 




'cutoff 
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«(r) = U a (r)V:(r) + CUt ° dE a U Ea (r)V£ a (r) . (1.9) 

0<E a <-\ J X 

As we have already mentioned in the introduction, most of the prac- 
tical calculations are done by solving Eq.(l.l) with a box boundary con- 
dition [2, 3, 4], i.e., by requiring that the solutions U a (r),V a (r) have a 
node at r = R (they are taken to be identically zero beyond R). This 
condition makes the spectrum of E a entirely discrete and allows the use 
of Eqs. (1.3-1.4) instead of Eq.(1.9). However, it is not clear how accu- 
rately one can mock up continuum effects like single-particle resonance 
contributions by using this discretization procedure. In order to avoid 
such ambiguities one needs solutions with proper boundary conditions. 
These solutions are illustrated in the next section for the case of a simple 
model [5, 6]. 

4. HFB SOLUTIONS: A SCHEMATIC MODEL 

In what follows we discuss the solutions of HFB equations in coordi- 
nate space for a mean field given by a square well potential, of depth Vq 
and radius a, and a pairing field, A, constant inside the same radius a 
and zero outside. We suppose also that the Fermi level, A, is given. For 
such a system the radial HFB equations inside the potential well, i.e. 
for r < a, are given by: 

,1 d 2 1(1 + 1) 2 . r2 

(-^2 r — + a ) U U- 6v l3 = 0. 

(i^r-^ + = 0. (1.10) 

where a 2 = f?(A + E + U ), (3 2 = |f (A -E + U ), 5 2 = §£A and 

Uq = — (Vq + V so l.s). The above equations have the following physical 
solution for any value of the quasiparticle energy : 

uij = A + ji(k + r) + A-ji(k-r) , 

vij = A + g + ji(k + r) + A_g_ji(k-r) , (1.11) 

where ji are spherical Bessel functions, k± = jfi{Uo + A ± \T(E 2 — A 2 )) 

and g ± = {E± ^(E 2 - A 2 ))/A. 

Outside the potential well both Uq and A are zero and the HFB 
equations are decoupled. In this case the type of solutions depends on 
the quasiparticle energy. Thus, for E < —A the solutions have the form: 

uij = Ah\ + \air) , 

Vlj = Bh\ + \(3 ir ) , (1.12) 
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where h\ are spherical Haenkel functions, a\ = %?(A + E) and 01 = 

^t(A — E). These solutions correspond to the bound quasiparticle spec- 
trum. 

For E > — A the spectrum is continuous and the solutions are: 

uij = C[cos(5i j )ji(air) - sin(<%)n/(air)] , 

vij = Dh\ +) (f3 ir ) , (1.13) 

where n\ are spherical Neumann functions and 5ij is the phase shift 
corresponding to the angular momentum (Ij). 

The constants entering in the wave functions above are fixed by the 
continuity conditions of the solutions and of their derivatives at r = a 
and the normalisation. In what follows we discuss only the continuous 
solutions, i.e., for energies E > —A. 

Of particular interest are the values of energies for which the wave 
functions have maximum localisations inside the potential well. These 
are the regions close to quasiparticle resonances, which can be defined 
as complex outgoing solutions of HFB equations. In HFB approach 
one distinguishes two types of quasiparticle resonances. One type cor- 
responds to the single-particle resonances of the mean field. For the 
pairing correlations an important role is played by those single-particle 
resonances which are close to the particle thereshold and have relatively 
high angular momentum. 

Another type of quasiparticle resonances corresponds to bound single- 
particle states. These resonant states, which appear due to the non- 
diagonal matrix elements of the pairing field, are specific of the HFB ap- 
proach. The resonant states corresponding to deep hole states have small 
widths. The states with very small widths can be eventually treated as 
quasibound states, normalized to unity in the same volume as the bound 
quasiparticle states. 

A special attention is paid usually to the continuum sl/2 states [6, 7]. 
Apart from the deep hole sl/2 states, which in HFB become narrow 
quasiparticle resonances with high quasiparticle energies, in drip line 
nuclei one may also find a loosely bound sl/2 single-particle state. In 
continuum HFB approach this state appears often as a broad quasi- 
particle resonance (see example below), close to the continuum edge, 
and its role in pairing correlations cannot be distinguished from the rest 
of non-resonant sl/2 continuum. This is different from a HF-BCS ap- 
proach, where the contribution of the loosely bound single-particle sl/2 
state to pairing correlations is not mixed with the rest of background 
single-particle continuum. 
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The structure of the continuum discussed above is essentially the same 
for general, self-consistent HFB calculations. In this case the HFB equa- 
tions for the continuum spectrum are integrated by starting far from 
the nucleus with the solution given by Eq.(1.13), which is propagated 
towards the matching point by the Numerov method. For each quasi- 
particle energy one calculates, by matching conditions, the phase shift. 
Then the energies (widths) of quasiparticle resonances are found from 
the energies where the derivative of the phase shift reaches its maximum 
(half of its maximum) value. This information is used afterwards to fix 
an appropriate energy grid (i.e., dense in the energy region of a reso- 
nance) for the calculation of the continuum contribution to the particle 
and pairing densities (see Eqs.(1.9)). 

5. HF-BCS APPROXIMATION 

The HF-BCS approximation is obtained by neglecting in the HFB 
equations the non-diagonal matrix elements of the pairing field. This 
means that in the HF-BCS limit one neglects the pairing correlations 
induced by the pairs formed in states which are not time-reversed part- 
ners. 

Particularly simple are the HF-BCS equations which include the effect 
of resonant continuum [8, 9]: 

A i = ^2(h i \ v \jJ) u j v j + ^2(i,i\V\ve l/ ,T7e^} / g°(e)u v (e)v v (e)(k , 

3 v " 

Me) = (9u(e)/9u(e))(J2(^ w \ V \jJ) u 3 v 3 

3 

+ ^2(iye l/ ,Ve^\V\i//e u ,,U7eV/) / gl,{el)u v i{€l)v v ,(el)del) 

ee (£(e)/g v (e))A v . (1.14) 

where g v {e) = 2 ^ +1 ^ff, 9v{ e ) ls the total level density and 5 V is the phase 
shift of angular momentum (l„j u ). In these equations the interaction 
matrix elements are calculated with the scattering wave functions at 
resonance energies and normalised inside the volume where the pairing 
interaction is active. The particle number condition is: 

= E«*+E / 9t(e)vl(e)de. (1.15) 

The energy factor g u (e) takes into account the variation of the localisa- 
tion of scattering states in the energy region of a resonance ( i.e., the 
width effect) and goes to a delta function in the limit of a very narrow 
width. For more details one can see Ref.[9]. 
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6. APPLICATION: CONTINUUM COUPLING 
IN DIFFERENT APPROXIMATIONS 

In order to illustrate the approximations discussed above we take as an 
example the nucleus 84 Ni, for which HFB calculations with discretized 
continuum can be found in the literature [4]. 



Table 1 Energies and widths of quasiparticle resonant states. In the two last columns 
the total occupancies for each (lj) channel are shown both in full and in box calcula- 
tions. 



j 


1 


E(MeV) 


r(MeV) 


Total occ. (full) 


Total occ. (box) 


1/2 





1.241 


0.69 






1/2 





20.956 


0.154 






1/2 





43.488 


10~ 6 


2.261 


2.299 


1/2 


1 


8.076 


0.36 






1/2 


1 


33.531 


0.101 


1.993 


1.994 


3/2 


1 


9.806 


0.60 






3/2 


1 


35.058 


0.104 


1.999 


2.016 


3/2 


2 


2.354 


0.641 






3/2 


2 


22.103 


0.066 


1.171 


1.183 


5/2 


2 


1.816 


0.076 






5/2 


2 


25.713 


0.006 


1.587 


1.529 


5/2 


3 


8.972 


0.919 


0.997 


1.006 


7/2 


3 


15.449 


2.070 


1.011 


1.015 


7/2 


4 


3.492 


0.027 


0.101 


0.110 


9/2 


4 


5.747 


0.001 


0.958 


0.952 


11/2 


5 


5.275 


0.049 


0.068 


0.079 



In all calculations discussed below the HF mean field is calculated 
with the Skyrme interaction SIII. The pairing force is given by Eq.(1.6) 
and the parameters are the ones used in Ref [4]. 

First, we present the results given by the continuum HFB equations 
solved with proper boundary conditions, as defined in the previous sec- 
tions. These results, referred below as "full" continuum calculations, 
are compared to the HFB calculations performed with box boundary 
conditions. 
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In Table 1 are shown the quasiparticle resonant states, of hole and 
particle type. For each (Ij) channel are shown also the total occupancy 
obtained in full and box calculations. One can notice the large widths of 
quasiparticle states corresponding to bound single-particle states close to 
the Fermi energy. The widths are obtained from the phaseshift behaviour 
around ir/2 except for the three states sl/2, pl/2 and p3/2 close to the 
Fermi energy for which the widths are extracted from the occupancy 
profiles. 

The results for the pairing field, pairing density and particle density 
are shown in Figs. (1-3). On can see that the box calculations generally 
overestimate the pairing correlations. This is also seen by comparing 
the Fermi levels (full: -1.026 MeV; box: -1.205 MeV), averaged gaps 
(full: 1.230 MeV; box: 1.448 MeV), pairing energies (full: -18.092 MeV; 
box: -24.081 MeV) and total binding energies (full: -651.079 MeV; box:- 
654.708 MeV). 

Next, we compare the results given by the HF-BCS approximation 
presented in the previous section, to the HFB calculations of Ref.[4]. In 
Ref. [4] the HFB equations are diagonalised in a basis formed by the 
single-particle states selected by a box of a finite radius. The single- 
particle continuum is cut at 5 MeV. This energy cutoff is different from 
the HFB box calculations presented above, where the energy cutoff is 
much higher. Up to 5 MeV, one finds three single-particle resonances, 
<i 3 / 2 , 57/2 an d h n / 2 , which are treated in the HF-BCS approximation 
given by Eqs.(1.14). The HFB and HF-BCS results for the pairing field 
are shown in Fig. 4. In the same figure is shown the HF-BCS result 
obtained by replacing each single-particle resonance by a discrete state, 
normalized in the same volume as the one used for the bound states. As 
seen from Fig. 4, this result is closer to HFB calculations. This similarity 
is due to the fact that in the HFB calculations each resonant state is 
actually represented by a unique state normalised in the box. When the 
effect of the width of resonant states is included in HF-BCS equations, 
the pairing correlations are decreasing. This trend is similar with the 
one observed in continuum HFB calculations presented above. Thus, 
one can conclude that, in order to take fully into account the effect of 
resonant continuum upon pairing correlations, one needs to solve the 
HF-BCS and HFB equations with proper boundary conditions. 
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Figure 1 Pairing fields in full (full line) and box (dashed line) calculations. 
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Figure 2 Pairing densities in full (full line) and box (dashed line) calculations. 
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Figure 4 Pairing fields calculated in different approximations: HFB of Ref. 
(stars), resonant HF-BCS with widths (solid) and without widths (dashed). 
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7. SUMMARY 

In these lectures, we have discussed how the continuum coupling and 
pairing correlations are calculated in drip line nuclei. The discussion was 
restricted to the HFB and HF-BCS approximations. In the first part, we 
have shown how one can construct in the HFB approach solutions with 
proper boundary conditions for the continuum spectrum. Then, for the 
case of a simple model, we have analysed the structure of quasiparticle 
continuum. In particular we have discussed the treatment of resonant 
continuum both in HFB and HF-BCS approaches. 

In the second part we have shown, for the case of a neutron rich nu- 
cleus, how different treatments of continuum can affect the pairing cor- 
relations. One concludes that the solutions based on discretized contin- 
uum (box boundary conditions) can overestimate the pairing correlations 
when compared with proper continuum HFB or HF-BCS calculations. 
This may have consequences on the predictions of drip lines. 
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